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ONR  GRANT  #N0001 4-06-1-1 369 
PI:  Roger  M.  Samelson 

Title:  Predictability  and  Dynamics  of  Geophysical  Fluid  Flows  -  GRA  Extension 

Research  under  this  grant,  an  extension  of  Grant  #NOOO  14-98- 1-08 13  that  supported  the 
completion  of  the  Ph.D.  thesis  of  Christopher  L.  Wolfe,  focused  primarily  on 
computations  of  unstable  nonlinear  periodic  solutions,  time-dependent  normal  modes 
(Floquet  vectors),  and  singular  vectors  in  a  two-layer  quasi-geostrophic  channel  model. 
The  model  was  studied  in  a  strongly  nonlinear  regime,  in  which  small  disturbances  to  an 
unstable,  steady,  zonal,  baroclinic  shear  flow  grow  to  finite  amplitude  and  continue  to 
vacillate  irregularly  for  arbitrarily  long  times.  The  computation  of  time-dependent, 
normal-mode  disturbances  to  unstable,  nonlinear,  time-periodic  basic  flows  in  a  high¬ 
dimensional  geophysical  fluid  model  opens  a  new  perspective  on  the  analysis  of 
disturbance  growth  in  time-dependent  flows,  and  on  the  closely  related  problem  of  error 
growth  in  predictive  models  of  time-dependent  flows. 

The  results  of  this  one-year  extension  confirm  and  extend  the  results  reported  for  the 
previous  Grant  #N00014-98-l-0813.  In  the  strongly  nonlinear  regime,  unstable  periodic 
solutions  were  identified  and  computed  using  an  efficient  Newton-Picard  method.  The 
complete  Floquet  and  singular  vector  spectra  were  computed  for  strongly  nonlinear 
cycles  and  the  dynamical  structure  of  the  resulting  modes  were  analyzed.  The  Floquet 
vectors  fall  into  two  classes  with  direct  physical  interpretations:  wave  dynamical  (WD) 
modes  and  damped-advective  (DA)  modes.  Singular  vectors  were  computed  and 
analyzed  in  terms  of  the  two  classes  of  normal-mode  solutions,  and  showed  a  similar 
dynamical  splitting.  A  new,  efficient  method  was  developed  for  the  recovery  of 
Lyapunov  vectors  from  singular  vectors,  which  promises  the  possibility  of  computing 
Lyapunov  vectors  directly  for  the  first  time  for  very  large  systems  such  as  numerical 
weather  prediction  models. 
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